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Abstract Current understanding of core collapse and thermonuclear supernovae is re- 
viewed. Recent progress in unveiling the nature of cosmic gamma-ray bursts 
(GRB) is discussed, with the focus on the apparent link of several GRBs with 
an energetic subclass of stellar explosions, type Ib/c core-collapse supernovae. 
This relation provides the strong case that the GRB phenomenon is connected 
with the final stages of massive star evolution and possibly with the formation of 
neutron stars and black holes. 
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Introduction 

Seventy years ago W. Baade and F. Zwicky (1934) (SJ were the first to 
point out tliat one of tlie briglitest astronomical phenomenon, supernovae stars 
(SNe), can be due to explosions of massive stars at the end of their evolution. 
The formation of a dense neutron core (neutron star) results in a sudden en- 
ergy release of order of the gravitational binding energy of the neutron star 
which amounts to Eg ~ —GMffg/Rj\fs ~ 10^^ ergs for the canonical val- 
ues of the NS mass M^vs ^ and radius R^s ~ 10 km. It was soon 
recognized by Gamow and Schoenberg (1941) [491 that most of this energy 
comes into neutrino emission. Twenty years after, Hoyle and Fowler (1960) 
showed that energy released in type la supernovae are connected to ther- 
monuclear burning of a degenerate stellar core. Here the available energy is 
~ 0.007M f« 2 X 10^^ ergs for the Chandrasekhar mass of the white dwarf. 

In the end of 1960s cosmic gamma-ray bursts (GRBs) were discovered by 
gamma-ray satellites ||79l|98l. Largely due to the inability of precise local- 
ization of the GRB position on the sky using gamma-ray detectors only, the 
origin of GRBs was as enigmatic as that of SNe before Baade and Zwicky's 
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suggestion until late 1990s, when the first successful localizations of GRBs us- 
ing their afterglow emission were made by BeppoSAX satellite in X-rays Q5||. 
The detection of X-ray afterglows several hours after GRB allowed dedicated 
follow-up observations of the GRB error boxes to be carried out using powerful 
optical (e.g. 1 1 18 1) and radio telescopes |40 1, in which rapidly decaying after- 
glow emissions were also detected. Quite soon after that, in the spring 1998, 
a bright peculiar nearby supernova 1998bw was found within the error box of 
GRB 980425 L45, .80J, suggesting the link between GRBs and SNe. Presently 
we have several unequivocal associations of GRBs with very energetic type Ibc 
supernovae called "hypernovae" (see below, Section|31l. 

Being connected to the evolution of stars, SN studies overlap with prac- 
tically all fields of the modern astronomy, from physics of tiny interstellar 
medium (e.g. to the formation of superdense neutron stars JHH- As was 
also pioneered by Baade and Zwicky, they are sources of astrophysical shocks 
in which cosmic ray particles are accelerated 1 17 1. 

Here we focus on some recent highlights in both core collapse and ther- 
monuclear supernova studies, which became possible mainly due to increas- 
ingly accurate radiation hydrodynamic calculations with a detailed treatment 
of neutrino processes. We also briefly describe recent success of asymmetric 
SN simulations (2D magneto-rotational collapse). Next we focus on recently 
established link between GRB explosions and energetic type Ibc supernovae 
(hypernovae) and discuss recent ideas on the GRB progenitors. We hypothe- 
size that different core collapse outcomes may lead to the formation of different 
classes of GRBs. 

1. Core collapse supernovae 

An extensive discussion of basic physics of core collapse supernovae can 
be found in fl2!i: the evolution of massive stars, core collapse, formation of 
stellar remnants and supernova nucleosynthesis are reviewed in II163I ; a recent 
concise discussion of problems and prospects for core collapse supernovae can 
be found in 1 1051 . 

In the end of thermonuclear evolution, the core of a massive star can lose 
mechanical stability for various reasons. In the stellar mass range 8Mq < 
M < 20Mq a partially degenerate core with mass close to the Chandrasekhar 
limit Mcore ~ ^ch and high density (p ~ 10^ — 10^" glow?) appears. Under 
these physical conditions, the chemical potential of degenerate electrons be- 
comes so high that neutronisation reactions + {A, Z) —>■ {A, Z — 1) + 
become effective even at zero temperature. At densities > 10^ g/cm^ de- 
generate electrons becomes relativistic so the adiabatic index of matter 7 = 
d\og P/dlog /9 — > 4/3, the critical value for loss of mechanical stability. Neu- 
tronisation of matter means its deleptonisation (decrease in the lepton number 
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Ye = Ne/Nij), SO the pressure at some moment increases slower than p^/^ 
though 7 is formally above 4/3 calculated at constant Yf,, and a catastrophic 
collapse begins. Temperatures are higher at larger stellar masses M > 20Mq 
so the collapse is initiated by photodissociation of nuclei (here 7 really be- 
comes < 4/3). For most massive stars with M > 60Mq pair creation makes 
7 < 4/3. 

The collapse occurs on the dynamical (free-fall) time scale tff ~ 1/ ^/Gp ~ 
a fraction of second. Adiabaticity holds so the entropy per baryon s = S/ks ^ 
const ~ 1 and even can increase due to non-equilibrium beta-processes. Low 
specific entropy (compared to H-buming phase at the main sequence, s 10 — 
15) prevents dissociation of nuclei until they "touch" each other at densities of 
the order of the nuclear density, /9„ ~ 2 x 10^"^ g/cm^. The collapse stops (if 
the core mass is below some Mmax), and bounce of the shock occurs at ~ 50 
km from the center. 

The bounce shock heats up deleptonized matter and rapidly spends most 
of its kinetic energy to destroy nuclei and produce plenty of free nucleons (n, 
p). Modified URCA-processes l38l becomes important: + p — > n + i^e> 

+n p+i'e and pair-neutrino annihilation takes place: e~^+e~ — > fe + z^e- 
At the typical collapse temperatures T ~ 10 MeV a lot of z^'s is produced 11681 
However, at densities p ~ 10^^ g/cm^ the mean free path of 10 MeV neutri- 
nos is by 5-6 orders of magnitude smaller than the size of the proto neutron 
star {R ~ 50 km) so the opaque "neutrinosphere" forms. Most of the core 
collapse neutrinos diffuse out of the neutrinosphere on a time scale ~ 10 sec- 
onds. First calculations of v spectra in core collapse SN were performed by 
D.K. Nadyozhin [108 109 | We should note that subsequent detailed calcula- 
tions (e.g. I lOlJ and references therein) did not change much these spectra. 
Thus, the modest ~ 10% fraction of the total neutrino energy released in the 
core collapse (~ 10^^ ergs) would be sufficient to unbind the overlying stellar 
envelope and produce the phenomenon of type II supernova explosion. 

Neutrino-driven explosions 

Thermal SN explosion mechanism was proposed by Colgate and White 1 32 1 . 
In this picture, part of the neutrino flux liberated in the core collapse is de- 
posited to the stellar mantle to make it unbound (~ 10^^ ergs is needed). Spe- 
cific mechanisms include neutrino-driven fluid instabilities, for example con- 
vection both above neutrinosphere and inside the proto-NS. Neutrino-driven 
convection, however, may not be as important as thought before, as follows 
from recent detailed 2D studies of convection |24|. Instead, other fluid insta- 
bilities such as newly found double-diffusive instability (the so-called "lepto- 
entropy fingers") l23ll . may effectively increase neutrino luminosity to help 
successful explosion. Note here that process + Ue ^ e'^ + ^7 + 7 
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may also be important above i^-sphere 0. This process was proposed as the 
energy source for GRB fireballs. 

It is now realized that even detailed ID-calculations of the core collapse 
supemovae fail to produce explosion (see e.g. reviews by f87l ll05l and refer- 
ences therein). The main reasons are that the bounce shock rapidly stalls (turns 
into a pure accretion front) over neutrinosphere at Rs ~ 200 km because of 
nuclear dissociations, and net neutrino heating of the freshly accreted material 
immediately after the front is insufficient. Burrows L25J notes that ID-models 
are about to produces a successful explosion, only 25% to 50% increase in the 
energy deposition rate inside the "gain region" is required. The stalled shock 
is not revived because the fresh fall-back matter is advected more rapidly than 
heated up by neutrinos. An additional heating by the viscous dissipation in 
differentially rotating models is considered in 11421 . Differential rotation of 
collapsing core seems quite plausible (see evolutionary calculations of rotat- 
ing stars |54 .55J). The energy stored in shear can amount to ~ 10^^ ergs 
for reasonably short 6 ms spin periods of proto neutron stars. The magneto- 
rotational instability (MRI) which arises in differentially rotating magnetized 
fluids with negative angular velocity gradients dQ/dr < 11491 l27l . MRI is 
thought to be responsible for turbulence in accretion disks and may operate in 
stars as well. This viscosity was found 11421 to be comparable with neutrino- 
driven convective viscosity and much larger than the neutrino viscosity. The 
account for additional energy dissipation due to the MRI viscosity in differen- 
tially rotating stellar cores allowed I142II to obtain a successful explosion in 
their ID-calculations. 

Transition from ID to multi-dimensional hydrodynamic calculations has en- 
abled a more accurate treatment of fluid instabilities. Presently, it is recognized 
that even with detailed account of Boltzmann neutrino transport, additional 
physical properties of the collapsing core, e.g. rotation, need to be coupled 
with multi-dimensional calculations to produce successful neutrino-driven ex- 
plosions (see e.g. f72lll05l ). 

Asymmetric explosions 

There are increasing observational evidence that supernova explosions are 
generically asymmetric. The indications are as follows. 

Spectroscopic observations. The asphericity of the explosion is directly 
inferred from spectropolarimetric observations |86|. The smaller the hydrogen 
mass of the presupernova envelope, the higher the degree of linear polarization 
is observed. In the canonical type Il-plateau supernovae with heavy hydrogen 
envelopes (SN1999em, 2003gd etc.) polarization is small (less than 0.5%), 
indicating low degree of asphericity. In contrast, in the stripped-envelope 
progenitors giving rise to type lib (SN1993J), type Iln (SN1998S), type Ic 
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(SN2003gf) or peculiar Ic (SN2002ap, maybe hypernova, see below) super- 
novae, th polarization degree is large, of order of 1-2%. The degree of polar- 
ization tends to increase with time (i.e. as more and more deeper layers of the 
ejecta become transparent). This suggests the relation of asymmetric super- 
nova explosion to pulsar kicks and, in extreme cases, to GRB explosions (see 
below). 

Pulsar kick velocities. Based on measurements of proper motions of ~ 
100 radio pulsars, it is well established now 1^ |7| that young neutron stars 
demonstrate space velocities up to several 100 km/s which cannot be explained 
without introducing additional natal kicks. 

Shape of young galactic SNR. The asphericity of SNRs 1987a, Cas A, 
N132D etc. is directly seen in multiwavelength observations ril6''6Tl. Recent 
Chandra X-ray observations of galactic SNR W49B [77 1 suggest this SNR 
may be the remnant of an asymmetric hypernova explosion in our Galaxy. 

Strong mixing of Ni-56. Strong mixing of Ni — 56 in the outer layers of the 
ejecta is required to explain the observed SN light curves (especially SN1987a, 

GoTl). 

Different mechanisms for the SN asymmetry have been proposed. 

Hydrodynamic mechanism. Large-scale neutrino-driven overturn between 
the shock and neutrinosphere can cause the SN anisotropy explosion |26 1. Pul- 
sar kicks up to 500 km/s from neutrino convection were obtained in recent 2D 
hydrodynamic calculations with Boltzmann neutrino transport |24|, while in 
3D-calculations (but with less rigorous treatment of neutrino transport) 1 43 1 
high initial velocities of proto-NS were found to be damped by neutrino emis- 
sion generally to less than 200 km/s. However, the postbounce stellar core 
flow is found (both in 2D and 3D) to be subjected to generic non-spherical per- 
turbations of the accretion shock leading to the development of the so-called 
"stationary accretion shock instability" l22l . This large-scale (/ = 1, 2 modes) 
fluid instability was also found in multi-dimensional simulations in Ref. |73|, 
and may be responsible for observed bimodal pulsar velocity distribution |7 1. It 
is important that the developed accretion-shock instability can result in bipolar 
explosions even in the absence of rotation. 

Magnetorotational mechanism. Rotation of the pre-supemova core and 
magnetic field increase during the collapse, and was proposed in the beginning 
of 1970s lfT3l as an alternative mechanism for core-collapse SNe. The magne- 
torotational mechanism by its nature is asymmetric and can launch antiparallel 
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jets during the explosion. First simulations were made in Ref. lEllIilD. Re- 
sults of recent 2D MHD calculations |3|51 are very encouraging: strong dif- 
ferential rotation in the presence of (maybe initially weak) magnetic field was 
shown to increase the magnetic pressure and form a MHD shock. As a result, 
rotational energy of neutron star is converted to the energy of the radial expan- 
sion of the envelope. The exponential growth of magnetic field in differentially 
rotating collapsing stellar cores due to MRI 1 1 1 appears to help the magnetoro- 
tational explosion as well. The magnetorotational explosion was found in [5^ 
to be essentially divided into three stages: linear growth of the toroidal com- 
ponent of the magnetic field due to twisting of the magnetic filed lines, the 
exponential growth of both toroidal and poloidal field components due to de- 
velopment of MHD instabilities, and the formation of a MHD shock leading to 
the explosion. These 2D MHD simulations have shown that the obtained 
energy of the magnetorotational explosion ~ 0.6 x 10^^ ergs is sufficient to 
explain type II and type lb core-collapse supernovae. 

Formation and subsequent collision of a close binary NS+NS system inside 
the stellar interiors. V.S. Imshennik 1 64 1 suggested the following scenario: 
a rapidly rotating core collapse results in the core fusion due to rotational insta- 
bility (the original idea goes back to von Weizsaecker |158|), then the binary 
NS system coalesces due to gravitational radiation losses. In this scenario, the 
lighter NS (with larger radius) first fills its Roche lobe and start losing mass 
as the orbit shrinks due to gravitational radiation. The NS gets unstable when 
M < Mmin ~ O.IMq providing an energy release of ~ lO^'^ erg in neutri- 
nos llT8ll33l l34l. This mechanism has been further elaborated in Ref. ll65l l2l. 
It was suggested for SN1987a explosion [66 1 to explain double LSD neutrino 
signal separated by ~ 5 hours observed from SN 1987a. A point of concern 
here may be a rapidly rotating pre-supernova core required for this mechanism 
to operate. Recent calculations by the Geneva group of stellar evolution with 
rotation but without magnetic field ll57l indicate that there are enough angular 
momentum in the pre-supernova core, while calculations with (even approxi- 
mate) account for magnetic field 1 55 1 show a severe (30-50 times) reducing of 
the final rotation of the collapsing iron core in massive stars. The core frag- 
mentation during collapse with rotation was also not found in 3D calculations 
in Ref. M . 

Non-standard neutrino physics. To explain high pulsar kicks, neutrino 
asymmetries in high magnetic fields pertinent to young pulsars have been in- 
voked [ 29l [141 183 1. The neutrino asymmetry is a natural consequence of the 
asymmetry of basic weak interactions in the presence of a strong magnetic 
field. Neutrino oscillations requiring a sterile neutrino with mass 2-20 keV 
and a small mixing with active neutrinos was proposed in Ref. ISTl . However, 
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a critical analysis in Ref. 1731 suggests that neither rapid rotation nor strong 
magnetic field of the young neutron star can be unequivocally inferred from 
observations of core-collapse SNR 1987 A and Cas A, so the solution to the 
problem of high pulsar kicks probably should be looked for within the frame 
of the conventional hydrodynamic SN explosion mechanism (see above). 

2. Thermonuclear supernovae 

Thermonuclear supernovae constitute a separate very important class of stel- 
lar explosions. In contrast to core-collapse SNe they demonstrate very similar 
light curves and which allowed their using as "standard candles" in modem 
cosmology (e.g. [128] and references therein). SN la are due to thermonuclear 
burning of a (C-i-0) white dwarf with M ~ Mch 1 59 1. A recent deep review of 
the SN la explosion models can be found in Ref. |56 1. In the modern picture 
(e.g. 11651 '). thermal instabilities in degenerate matter with p ~ 2 — 9 x 10^ 
g/cm^, T ~ 7 X 10^ K, after ~ 1000 years of core convection initiate flame 
ignition within a typical ignition radius ~ 150 — 200 km and could result in 
the complete destruction of the white dwarf. In principle, at higher density a 
neutron star could be formed due to electron-capture processes (via accretion- 
induced collapse of a 0-i-N-i-Mg core or accretion onto a 0-Ne-Mg white dwarf 
in a binary system 1 133. .621 ). but higher densities are disfavored from evolu- 
tionary viewpoint. 

Possible scenarios for the formation of type la SNe include double degener- 
ate white dwarf mergings or accretion onto a massive white dwarf in a symbi- 
otic binary (see discussion of different formation channels for SN la progeni- 
tors in Refs. Ill66lll3l1ll67l ). Sub-Chandrasekhar mass models for SNla due 
to an external trigger (e.g. detonation of the accreted He layer) may result in 
subluminous SN la (like 1991bg). A SN 1.5 model in which the degenerate 
core explodes at the late asymptotic giant phase of the evolution of an inter- 
mediate mass star was suggested in Ref. II63I SN 2002ic and SN 1997cy in 
whose spectra hydrogen absorption lines were discovered probably belong to 
this class ll30l . In this model the hydrogen envelope ejection is synchronized 
(within ~ 600 years) with the explosion of a contracting white dwarf as its 
mass approaches the Chandrasekhar limit. 

After the initial thermonuclear ignition of the degenerate core interior, a 
strong temperature dependence of the nuclear reaction rates (oc T^^ at ^ MeV 
temperatures) leads to the formation of very thin burning layers propagating 
conductively with subsonic speeds (deflagration, i.e. a flame) or burning due 
to shock compression (supersonic detonation). Both propagation modes are 
linearly unstable to spatial perturbations and presently are treated using multi- 
dimensional calculations. The prompt detonation of the degenerate core as SN 
la mechanism was first studied by D. Arnett |6|. However, the mechanism 
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is inconsistent with observations as too little intermediate-mass elements are 
bom in detonation, and was found not to operate because the core at ignition 
is insufficiently isothermal |T60|. Deflagration to detonation (delayed detona- 
tion) burning regime was suggested in Ref . 1 67 1 and further elaborated in 1 78 1 . 
Pure carbon deflagration with convective heat transport was proposed in Ref. 

nioi . 

The main problem encountered in studies of the SN la explosion mecha- 
nisms can be formulated as follows: the prompt detonation produces enough 
energy for explosion (thermonuclear runway energy a white dwarf ~ 0.007Mc^ ~ 
1.5 X 10^^ ergs) but gives incorrect nucleosynthesis yields (mainly iron peak 
elements, in contrast with observations showing significant abundance of the 
intermediate mass elements), the pure deflagration is too slow and must be ac- 
celerated. The acceleration of the deflagration front is achieved by involving 
different flame instabilities: Landau-Darrieus flame instability, which leads 
to the flame front fractalization - wrinkles and folds, so the surface area of 
the flame effectively increases (see analytical treatment in fT9ll and numerical 
simulations in |129|); (2) Rayleigh-Taylor (RT) instability resulting from the 
buoyancy of hot, burned material in the dense, unburned surroundings. The 
RT instability generically develops after the initial ignition in degenerate cores 
turbulent deflagration and leads to the appearance of the hot bubbles ("mush- 
rooms") floating upwards while spikes of cold fluid falls down. Secondary 
hydrodynamic instabilities due to shear along the bubble surface rapidly leads 
to the turbulence. The key role of the RT instability in the presently concurrent 
turbulent deflagration and delayed-detonation models for SN la explosions is 
fully confirmed in 3D numerical simulations fTzT.Wf]. 

Recent progress in the physics of the SN la explosions has been done mainly 
due to multi-dimensional numerical simulations of the flame propagation in the 
degenerate star. For example, multidimensional Chandrasekhar mass deflagra- 
tion simulations [127| indicate acceleration of the turbulent combustion front 
up to 30% of the speed of sound, which is enough to produce an explosion 
without transition to detonation. Generally, 3D-models are found to be more 
energetic. Multi-dimensional calculations of nucleosynthesis in the pure defla- 
gration Chandrasekhar model |144J shows that turbulent flame converts about 
50% of carbon and oxygen to ash with different composition depending on the 
density of the unburned material. To burn most of the material in the center (so 
that to avoid unobserved low-velocity carbon, oxygen and intermediate-mass 
elements in the spectra), the model requires a large number of ignition spots. 
In contrast, 3D simulations of the delayed-detonation SN la explosion BHl in- 
dicate that there is no such problem in this model and the explosion energy is 
higher than that obtained in the pure deflagration burning. It remains to be seen 
from future observations which model is more adequate. 
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Note that uncertainties in the ignition conditions of the degenerate star leads 
to some irreducible diversity of the explosion kinetic energy, peak luminosity, 
nickel production for the same initial configuration. Modeling of light curves 
of SN la turns out to be a powerful tool to check the SN la explosion models 
(see recent calculations by multi-group radiation hydrocode STELLA (ITl ). 

3. Gamma-ray bursts 

GRBs have remained in the focus of modem astrophysical studies for more 
than 30 years. After the discovery of GRB afterglows in 1997 |35 1, the model 
of GRB as being due to a strong explosion with isotropic energy release of 10^^ 
ergs in the interstellar medium became widely recognized. Various aspects 
of GRB phenomenology are discussed in many reviews: observational and 
theoretical studies are summarized in 1,60.1 . first observations of afterglows are 
specially reviewed in I119L GRB theory is extensively discussed in I1021ll69ll . 

A widely used paradigm for GRBs is the so-called fireball model (e.g. re- 
views by Piran 1 120, 121 1 and references therein). In this model, the energy 
is released in the form of thermal energy (its initial form is usually not spec- 
ified) near the compact central source (at distances and is mostly converted 
into leptons and photons (the fireball itself). The relativistic outflow (wind) is 
formed driven by the high photon-lepton pressure (generically in the form of 
two oppositely directed narrow collimated jets) I113II1351 . The fireball inter- 
nal energy is converted to the bulk motion of ions so that relativistic speed with 
high Lorentz-factors (typically, F > 100) is achieved during the initial stage 
of the expansion; the ultrarelativistic motion is in fact dictated by the need to 
solve the fireball compactness problem (see l20ll for a detailed discussion and 
references). The kinetic motion of ions is reconverted back into heat in strong 
coUisionless relativistic shocks at typical distances of 10^^ cm. Assuming the 
appropriate turbulence magnetic field generation and particle acceleration in 
the shocks, energy thermalized in the shocks is emitted via synchrotron and 
inverse-Compton radiation of shock-accelerated electrons 1 125 1 (see 1 156| for 
a review), which is identified with the GRB emission. A shell of ultrarelativis- 
tically moving cold protons produces a blast wave in the surrounding medium, 
forming an external shock propagating outward and reverse shock that propa- 
gates inward and decelerates the explosion debris. Most energy of the explo- 
sion is now carried by the external shock which decelerates in the surrounding 
medium. Assuming magnetic field generation and particle acceleration in the 
external shock, the afterglow synchrotron emission of GRB is produced in ra- 
dio fTBl . optical |75 1041 and X-rays |151, 104 1. Note that at this stage 
the memory of the initial explosion conditions is cleaned, and the dynamical 
evolution of the external shock is well described by the Blandford-McKee self- 
similar solution 1 15 1, a relativistic analog of the Sedov-von-Neumann-Taylor 
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solution for strong point-like explosion. This explains the apparent success in 
modeling the GRB afterglow spectral and temporal behavior in the framework 
of the simple synchrotron model II159I . irrespective of the actual nature of the 
GRB explosion. 

There is no consensus thus far about the origin of the GRB emission itself. 
Within the fireball model, the GRB can be produced by internal fireball dissi- 
pation (the internal shock wave model, e.g. [ 126]), or in the external blast wave 
decelerating in the ambient (inhomogeneous) medium 1.103.. 37,1 . The fireball 
model is known to face some important problems (for example, baryon con- 
tamination of the fireball, the microphysics of magnetic field generation and 
particle acceleration in collisionless ultrarelativistic shocks etc., see a critical 
review in |93 36 1). In Ref. [93 1 an alternative to the fireball model is analyzed 
in which large-scale magnetic fields are dynamically important. Whether the 
GRB jets are hot (fireball model) or cold (electromagnetic model) remains to 
be determined from future observations. Here crucial may be spotting the very 
early GRB afterglows and measuring polarization of prompt GRB emission 
(see |92| for the short-list of the electromagnetic model predictions). Note 
that irrespective of the mechanism mediating the energy transfer from the cen- 
tral source to the baryon-free region, many essential features of the observed 
non-thermal GRB spectra can be reproduced in some general physical models 
of prompt gamma-ray emission of GRBs, e.g. by synchrotron self-Compton 
emission of plasma with continuously heated nearly monoenergetic electrons 

An important open issue is whether GRBs can be the sources of ultra- 
high energy cosmic rays (UHECR). This association was suggested in Refs. 
1 155 106, 150| based on similarity of the energy release in GRBs ~ 10^^ erg 
per year per cubic Megaparsec with what is observed in UHECRs and assum- 
ing effective proton acceleration in relativistic collisionless shocks. The mech- 
anism of UHECR production in GRBs is still uncertain, but basic requirements 
for proton acceleration to high energies in mildly -relativistic shocks (pertinent 
to the internal shock model of GRBs) appear to be satisfied. See Ref. 1 1571 for 
more detail and discussion. See also lectures by M. Teshima and M. Ostrowski 
on this School. 

Below we focus on the observed association of GRBs with an energetic sub- 
class of core-collapse supernovae, type Ibc SNe, which with each new finding 
provides an increasing evidence that the GRB phenomenon is related to the 
evolution of most massive stars and formation of relativistic compact objects 
(neutron stars and black holes). 
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4. Supernova - GRB connection 
Theoretical grounds: the collapsar model 

The connection of GRBs with stellar explosions was first proposed theoret- 
ically. Woosley (1993) [161] considered a model of accretion onto a newly 
formed rotating black hole to power the GRB fireball. The progenitor to GRB 
in this model is a rapidly rotating Wolf-Rayet (WR) star deprived of its hydro- 
gen and even helium envelop due to powerful stellar wind or mass transfer 
in a binary system. Dubbed by Woosley himself as "failed type lb super- 
novae", this model is now called the collapsar model f94l. In this model, a 
massive (> IbM^) rotating star with a helium core > lOM© collapses to 
form a rapidly rotating BH with mass > 2 — 3Mq. The accretion disk from 
the presupemova debris around the BH is assumed to be the energy source for 
GRB and is shown to be capable of providing the prerequisite 10^^ — 10^^ ergs 
via viscous dissipation into neutrino-antineutrino fireball. The energy released 
is assumed to be canahzed in two thin antiparallel jets penetrating the stellar 
envelope. 

Another possible energy source in the collapsar model could be the electro- 
magnetic (Poynting-dominated) beamed outflow created via MHD processes, 
much alike what happens in the active galactic nuclei powered by accretion 
onto a supermassive BH. The estimates show that the Blandford-Znajek (BZ) 
process 1 16 1 in the collapsar model (e.g. 1 85 1) can be a viable candidate for the 
central engine mechanism for GRBs, provided somewhat extreme values for 
BH spin (the Kerr parameter a = Jc/GM"^ ~ 1) and magnetic field strength 
in the inner accretion disk around the BH {B ~ 10^^ — 10^^ G). In that case 
the rotating energy of BH (up to 0.29Mb/iC^ for a = 1) is transformed to the 
Poynting-dominated jet with energy sufficient to subsequently produce a GRB. 

Yet another source of energy in the collapsar model could be the rotation 
energy of a rapidly spinning neutron star with high magnetic field (magnetar), 
as originally proposed by Usov [146 147 1. As in the BZ-based models, the 
GRB jets are Poynting-dominated. Lyutikov and Blandford l93l develop the 
electromagnetic model, which postulates that the rotating energy of the GRB 
central engine is transformed into the electromagnetic energy (for example, 
in a way similar to the Goldreich-Iulian pulsar model) and is stored in a thin 
electromagnetically-dominated "bubble" inside the star. The bubble expands 
most rapidly along the rotational axis, breaks out of the stellar- envelopes and 
drives the ultrarelativistic shock in the circumstellar material. In contrast to the 
synchrotron GRB model, here GRB is produced directly by the magnetic field 
dissipation due to current-driven instabilities in this shell after the breakout. 
The energy transfer to GRB is mediated all the way by electromagnetic field 
and not by the ion bulk kinetic energy. It remains to be checked by observations 
whether the EM or fireball model for GRB emission is correct. 
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A different scenario is the so-called "supranova" model for GRB proposed 
in 11521 involves a delayed collapse to BH through the formation of a super- 
massive rotationally supported neutron star (see also [153] for another vari- 
ant of this model). In the supranova scenario GRB is associated with the BH 
formation and occurs in baryon-clean surrounding after the initial supernova 
explosion. The original scenario predicts several weeks - month time delay 
between the SN explosion and GRB. The critical comparison of the supranova 
mechanism with observations (showing its inconsistency) can be found in |36|. 
In this paper an interesting extension of the supranova scenario is considered 
in which the second collapse occurs minutes or hours after the primary SN 
explosion. In this model, a rotationally-supported magnetized neutron star is 
formed in a SN type Ic explosion and rapidly loses energy via powerful pulsar 
wind along the rotational axis. The wind drills the baryon-clean polar cones 
through which a newly formed rapidly spinning BH generates relativistic jets 
producing a GRB (for example, by BZ mechanism). Though no rapid rotation 
of magnetized neutron star has been obtained in stellar evolution calculations 
(nor can it definitely be inferred from the existing observations, cf. ll73l ). the 
spin increase with the remnant's mass found in Ref. [55 1 may provide credence 
to this mechanism. Some signatures from the pulsar-wind heated supernova 
shell can be tested in future observations. 

Observational evidence: GRB -supernova associations 

First hint on the association of GRBs with SNe came from the apparent time 
coincidence (to within about a day) of GRB 980425 with a peculiar supernova 
SN 1998bw |45|. SN 1998bw occurred in a spiral arm of nearby (redshift 
z = 0.0085, distance ~ 40 Mpc) spiral galaxy ESQ 184- G82. Such a close 
location of GRB 980425 rendered it a significant outliers by (isotropic) energy 
release AEiso ~ 10^*^ erg from the bulk of other GRBs with known energy 
release, and even from a beaming-corrected mean value of GRB energies of 
~ 10^^ erg |41|. 

Now the most convincing evidence for GRB-SN association is provided by 
spectroscopic observations of expanding photosphere features in late GRB af- 
terglows. Especially strong is the case of a bright GRB 030329 associated with 
SN 2003dh [58,138,^^, 76 J. Spectral observations of the optical afterglow 
of this GRB revealed the presence of thermal excess above non-thermal power- 
law continuum typical for GRB afterglows. Broad absorption troughs which 
became more and more pronounced as the afterglow faded indicated the pres- 
ence of high-velocity ejecta similar to those found in spectra of SN 1998bw. 
Despite these strong evidences, there are some facts which cannot be explained 
by simple combination of the typical SN Ibc spectrum and non-thermal power- 
law continuum. For example, the earliest spectroscopic observations of GRB 
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030329 of optical spectra taken on the 6-m telescope SAO RAS 10-12 hours af- 
ter the burst I136II showed the presence of broad spectral features which could 
not be produced by a SN at such an early stage. The complicated shape of the 
optical light curve of this GRB with many rebrightenings 1 88 1 and polarization 
observations made by VLT 1 52 1 suggest a clumpy circumburst medium and re- 
quire additional refreshening of shocks (if one applies the synchrotron model, 
e.g. |51|). 

Another interesting example of GRB-SN connection is provided by GRB 
031203. This GRB is one of the closest (z = 0.105) known GRBs and is found 
to be intrinsically faint, /\Eiso ~ 10^° ergs ( llT3ti ni^ ^ The low energy 
release in gamma-rays is confirmed by the afterglow calorimetry derived from 
the follow-up radio observations 1 137 1 and allows this GRB to be considered as 
an analog to GRB 980425. It is important that the low energy release in these 
bursts can not be ascribed to the off-axis observations of a "standard" GRB 
jet (unless one assumes a special broken power-law shape of GRB luminosity 
function, see ||53|). However, a bright type Ib/c supernova SN 20031w was 
associated with GRB 031203 as suggested by the rebrightening of the R light 
curve peaking 1 8 days after the burst and broad features in the optical spectra 
taken close to the maximum of the rebrightening [31 143 96 46|. 

The comparison of radio properties of 33 SNe type Ib/c with those of mea- 
sured radio GRB afterglows allowed Berger et al. fTO] to conclude that not 
more than few per cents of SNe type Ib/c could be associated with GRBs, 
which explains the observed small galactic rate of GRBs. However, it still re- 
mains to be studied how much intrinsically faint GRBs like 980425 and 03 1203 
can contribute to the total GRB rate. 

5. Hypernovae 

Core-collapse supemovae with kinetic energy of the ejecta ~ 10 — 30 times 
as high as the standard 1 foe (Ifoe = 10^^ erg) are now collectively called 
"hypernovae". The term was introduced by B. Paczynski shortly after the dis- 
covery of first GRB afterglows in 1997 by the Beppo-S AX satellite 1 1 14 1 based 
on qualitative analysis of possible evolutionary ways leading to cosmic GRB 
explosions. 

SN 1998bw was exceptionally bright compared to other Ib/c SNe (the peak 
bolometric luminosity of order 10^^ erg/s, comparable to the SN la peak lumi- 
nosities). This points to the presence of a substantial amount of ^^Ni isotope, 
the radioactive decay thereof being thought to power the early SN light curves. 
The spectra and light curve of SN 1998bw was modeled by the explosion of 
a bare C-i-0 of a very massive star that has lost its hydrogen and helium en- 
velopes with a kinetic energy more than ten times typical SNe energies li68 J . 
and they called SN 1998bw a hypernova. 
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Since then several other SNe were classified as SN 1998bw-like hyper- 
novae by their spectral features and light curves: SN 1997ef, SN 2002ap, SN 
2003dh/GRB030329, SN 20031w/031203. Recently, SN 1997dq was dubbed a 
hypernova by its similarity with SN 1997ef 1 100|. 

Extensive numerical modeling of light curves and spectra of hypemovae 
(see fTll I for a recent review) confirmed the need of atypically high for core- 
collapse SNe mass of nickel (~ 0.1 — O.SMq) to be present in the ejecta in 
order to explain the observed hypernova properties. The rapid rise in of the 
observed light curves of the "canonical" SN 1998bw requires a substantial 
amount of ^^Ni to be present near the surface. This strongly indicates the 
important role of mixing during the explosion as nickel is synthesized in deep 
layers during a spherical explosion. This fact can serve as an additional evi- 
dence for non-spherical type Ic explosions. As we already stressed above, the 
asphericity appears to be a ubiquitous feature of core-collapse supernovae in 
general, culminating in bipolar hypernova explosions associated with GRBs. 

Spectral modeling suggests I112II that the broad-band spectral features gen- 
erally seen in early and maximum light of hypernovae signal very rapid photo- 
spheric expansion. In particular, authors of Ref. 1 112| notice the very unusual 
for other SNe fact that 01 (A = 7774^1) and Call IR (at A ~ 8000A) absorp- 
tion lines merge into a single broad absorption in early spectra of SN 1998bw, 
which indicates a very large velocity of the ejecta (the line separation 30000 
km/s). 

In general, varying (a) the progenitor C-i-0 core mass from 2 to ^ 14 solar 
masses, choosing (2) the appropriate mass cut (corresponding to the mass of 
the compact remnant, a neutron star or black hole Mc = 1.2 — AMq), and (3) 
mass of ^^Ni isotope (r^ 0.1 — 0.5Mq) and its mixing allow 11121 to reproduce 
the observed spectra and light curves of hypemovae. 

The analysis of nucleosynthesis in hypernovae suggests a possible classifi- 
cation scheme of supernova explosions 1.1 ILL In this scheme, core collapse in 
stars with initial main sequence masses M^s < 25 — 30Mq leads to the for- 
mation of neutron stars, while more massive stars end up with the formation 
of black holes. Whether or not the collapse of such massive stars is associated 
with powerful hypernovae ("Hypernova branch") or faint supernovae ("Faint 
SN branch") can depend on additional ("hidden") physical parameters, such as 
the presupernova rotation, magnetic fields. [39], or the GRB progenitor being 
a massive binary system component 1 145 , 1 17 1. The need for other parameters 
determining the outcome of the core collapse also follows from the continuous 
distribution of C-i-0 cores of massive stars before the collapse, as inferred from 
observations, and strong discontinuity between masses of compact remnants 
(the "mass gap" between neutron stars and black holes) |28 p. 

The mass of ^^Ni synthesized in core collapse also appears to correlate with 
M„,. In ordinary SNe (like 1987a, 1993j, 19941), Mm = 0.08 ± O.OSMq, 
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but for hypernovae this mass increases up to ~ O.SM© for the most energetic 
events. Large amount of ^^Ni in hypernovae suggests a different nucleosyn- 
thesis event. It was shown that unlike conventional core-collapse shock nu- 
cleosynthesis |I1621, nucleosynthesis in bipolar supernova explosions |95| or 
in relativistic modest-entropy massive wind from accretion disk around a BH 
11241 can in principle give the observed high amount of nickel in hypernovae. 
Another important consequence of hypernova nucleosynthesis can be larger 
abundances (relative to the solar one) of Zn, Co, V and smaller abundance of 
Mn, Cr, the enhanced ratios of a-elements, and large ratio of Si, S relative to 
oxygen (see fill! for further detail). 

It is also necessary to note that energy requirements for hypernova explo- 
sions ( A^; > 2 X 10^^ ergs) can hardly be provided by the most elaborated 
delayed neutrino explosion mechanism. Indeed, the net explosion energy in 
this mechanism comes mostly from nuclear recombination of matter inside the 
gain radius. Analytical |70| and numerical calculations LZl I suggest this mass 
to fall within 0.01 — O.IM© range, so assuming the recombination energy re- 
lease 8 MeV per nucleon results in ~ 10^^ ergs of the explosion energy (see 
also calculations of the SN explosion energy as a function of the progenitor 
mass in Ref. |42|). The magnetorotational mechanism |^ fails to produce the 
hypernova energies either (see above). Unless some non-standard physics does 
operate, the hypernova energies can be recovered from accretion of the rotat- 
ing collapse debris onto BH (AE' ~ 0.06 - 0A2AMaC^ depending on the BH 
spin) or by BZ mechanism (AE < 0.29M;,/iC^). In both cases, rapid rotation 
(and extremely large magnetic field for BZ to operate) of the presupernova is 
required. Here realistic self-consistent calculations still have to be done, which 
is a very difficult task. 

6. Progenitors of GRBs 

The GRB-SN connection leads to the generally accepted concept that mas- 
sive stars that lost their envelopes are progenitors of long GRBs (this limitation 
is due to the fact that predominantly long GRBs with duration > 2 s can be well 
localized on the sky and provide rapid alerts for follow-up multiwavelength 
observations). For short single-pulsed GRBs (a quarter of all bursts, see e.g. 
catalog by Stem et al. 1 139 J ) the binary NS-i-NS/NS-i-BH merging hypothesis 
llT8lll30ll69l remains viable (see also recent general relativistic hydrodynamic 
models of the launch and propagation of relativistic jets due to thermal energy 
deposition near the center of binary mergers |3 1). 

As we already noted, the emerging empirical evidence is that there ex- 
ist intrinsically faint, single-pulsed, apparently spherically-symmetric GRBs 
(980425, 031203) associated with strong hypernovae. These hypernovae re- 
quire maximal amount of nickel to be synthesized in explosion and large ki- 
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netic energies. On the other hand, another unequivocal hypemova SN 2003dh, 
associated with the "classical" GRB 030329, can be modeled with excep- 
tionally high kinetic energy (4 x 10^^ ergs) but smaller amount of nickel 
(~ O.35M0) and smaller mass of the ejecta (8 — IOMq) |99|. These param- 
eters were obtained assuming spherical symmetry, which is of course not the 
case for GRB 030329^ . But if this tendency is real and will be confirmed by 
later observations, we can return to our hypothesis [123 1 that there should be 
distinct classes of GRBs according to what is the final outcome of collapse of 
the CO-core of a massive star. If the collapse ends up with the formation of a 
neutron star, an intrinsically faint smooth GRB could be produced and a heavy 
envelope is ejected in the associated SN Ib/c explosion. The GRB energy in 
this case can be essentially the rotation energy of the neutron star ~ 10^^ — lO^'^ 
ergs, as in the electromagnetic model 1 146 1. If a BH is formed, a lighter en- 
velope is ejected with accordingly smaller amount of nickel and possibly with 
higher kinetic energy of the ejecta, and more energetic, highly variable GRB 
with a "universal" jet structure 11221 emerges fed by non-stationary accretion 
onto the BH. 

The GRB energy can be also interpreted in more exotic way requiring a 
new physics. The possible relation of GRBs to mirror dark matter was dis- 
cussed in Ref. l20l . The conversion of light axions from SNe to photons as 
the source of the GRB fireballs was considered in ll89l (TTl . Recently it was 
suggested |50| that ultramassive axions in the mirror world with the Peccei- 
Quinn scale fa ~ 10^ — 10^ GeV and mass 1 MeV can be produced in 

the gravitational collapse or in the merging of two compact stars. The axions 
tap most of the released energy and can decay ~ 1000 km away mostly into 
visible electron-positron pairs (with 100% conversion efficiency) thus creating 
the primary GRB fireball. The estimates show that successful short GRBs can 
be obtained in compact binary coalescences, while long GRBs can be created 
in coUapsars. In extended SN II progenitors, this energy may help the mantle 
ejection. In compact CO-progenitors for SN Ib/c axions decay inside the star, 
so depending on the stellar radius weaker or stronger GRBs associated with 
SNe type Ib/c explosions can be observed. In this picture again the collapse 
with the formation of a neutron star or BH may have different signatures. 

7. Conclusions 

Cosmic explosions including various types of supernovae and GRBs are nat- 
ural consequence of stellar evolution. The big efforts of different groups to the- 
oretically understand the physical mechanism(s) of core collapse SNe appear 
to be approaching the final phase. The thermal (neutrino-driven) mechanism 
for core collapse SNe is mostly elaborated at present but still fails to produce 
a strong explosion. There is understanding why stars do not explode by this 
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mechanism, and the road map is designed how to obtain the explosion. Re- 
cent multidimensional simulations with accurate incorporation of Boltzmann 
neutrino transport indicate that effectively a modest boost in the neutrino lumi- 
nosity is still required. This can be due to inclusion of new fluid instabilities 
and more accurate treatment of neutrino processes and microphysics (the neu- 
tron star equation of state) into the full 3D hydrodynamic calculations. The 
rotation and magnetic fields is not yet fully included into calculations, which 
is a challenging task. First results of 2D MHD calculations of the magne- 
torotational supernova explosion seem to be encouraging. The most energetic 
(hypemova) explosions, however, require an additional to neutrino source of 
energy, and the rotation and magnetic fields can be the principal ingredients. 

An impressive progress has been done in multidimensional calculations of 
thermonuclear explosions of degenerate dwarfs for type la supernovae. It is 
still however unclear whether pure deflagration or delayed detonation is at 
work in SN la. The important problem is to more precisely determine the ini- 
tial ignition conditions. Detailed radiation hydrodynamic modeUng revealed 
that SN la light curves proved to be very sensitive to the explosion models and 
thus can be used to check the models. 

In contrast to SNe, the nature of cosmic GRBs remains unclear. The most 
important recent progress in understanding GRBs was in establishing the link 
of at least part of them to unusually energetic SN Ibc explosions (hypernovae). 
At present several unequivocal GRB-SN associations are known. The two 
closest GRBs discovered so far (GRB 989425 and GRB 031203) proved to 
be intrinsically weak compared to the bulk of other GRBs with measured red- 
shifts. They both show a single-peak smooth gamma-ray hght curve with no 
signs of jet-induced breaks in the afterglows. In the third (most strong) case 
of the GRB-SN association, GRB 030329/SN 2003dh, the GRB light curve 
is two-peaked, the afterglows show evidence for jet. Modeling of the under- 
laid hypernovae light curve and spectra revealed the first two cases to require 
smaller kinetic energies but higher mass of the ejecta and the amount of the 
synthesized nickel than SN 2003dh. We tentatively propose that the tendency 
"weaker, more spherically symmetric GRB - stronger hypemova" may indi- 
cate the formation of a NS in the case of weak GRBs and of a BH in the case 
of strong variable GRBs as the final outcome of the core collapse. In the NS 
case the GRB energy comes from the rotational energy of neutron star and is 
possibly mediated by the electromagnetic field. When BH is formed the GRB 
energy source is the gravitational energy released during non-stationary accre- 
tion onto the black hole or the black hole rotation. It is not still excluded that 
the GRB phenomenon signals some new physics underlying the formation of 
compact stars. 

Observations of various types of supernovae in other galaxies and especially 
of a (long awaited!) galactic event by all available means (including neutrino 
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and gravitational wave detectors) should undoubtedly be crucial for further 
understanding physical mechanisms of cosmic explosions. We are sure that the 
increasing statistics of GRB/SNe in the nearest future obtained with new GRB- 
dedicated space missions like SWIFT will tell us much more on the nature of 
GRBs and their progenitors. 
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Notes 

1. A bright soft X-ray flux was inferred from XMM observations of evolving X-ray halo for this burst 
11481 . making it an X-ray rich GRB i,j g4 j ; this point of view was questioned by 11321 . 

2. Recent radial velocity measurements of the companion star in low-mass X-ray binary 2S 0921-630 
limits the mass of the compact object within the range 1.9 it O.25M0 < Mx < 2.9 ± 0.4Mq, making 
it the plausible high-mass neutron star or low-mass BH (731. A sUghtly higher mass interval 2.OAf0 < 
Mx < 4.3Af0 (1-sigma) was obtained in fl34'. 

3. For example, two-component modeling of SN 2003dh as a slowly moving high-mass equatorial 
ejection and almost discontinuous low-mass polar outflow 1 164 requires the high mass of synthesized Ni- 
56 ~ O.5Af0, as in other hypemovae 
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